// Loesung_von_Aufgabe_3.7_4_zwei_Massen

// Äquipotenziallinien im Bereich zweier Massen

float V; // Gesamtpotenzial 
float m1 = 1.4E12; // Masse 1
float m2 = 6.0E11; // Masse 2
float f = 6.6726E-11; // Gravitationskonstante

void setup()
{
  size(400, 300);
  background(255);
}

void draw()
{
  translate(width/2, height/2);

  // Die beiden Massen werden gezeichnet
  noStroke();
  fill(0);
  ellipse(-55, 0, 30, 30);
  ellipse(70, 0, 15, 15);

  /* Das Gravitattionspotenzial im Bereich der beiden Massen wird berechnet 
   und die Äquipotenziallinien werden gezeichnet */
  for (int j = -150; j < 150; j++)
  {
    for (int i = -200; i < 200; i++)
    {
      V = -f*m1/sqrt(pow(i+55, 2) + pow(j, 2) + 0.001) - f*m2/sqrt(pow(i-70, 2) + pow(j, 2) + 0.001);

      if (V >= -1.0-0.05 && V <= -1.0+0.05 || V >= -1.5-0.05 && V <= -1.5+0.05 || V >= -2-0.05 && V <= -2+0.05 || V >= -3-0.05 && V <= -3+0.05)
      {
        stroke(255, 0, 0);
        strokeWeight(5);
        point(i, j);
      }
    }
  }
}